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ABSTRACT 

We investigate the implications for neutrino mixing implied by the results of all 
neutrino oscillation experiments and by the standard Big-Bang Nucleosynthesis 
constraint on the number of light neutrinos. 

Many experiments searching for neutrino oscillations have been done in the last 
30 years using solar, atmospheric, reactor and accelerator neutrinos. The majority 
of these experiments reported a negative result, but there are three positive indi- 
cations in favor of neutrino oscillations coming from the results of solar neutrino 
experimentJ, of atmospheric ueutrino experiment J arrd of the LSND accelerator 

~^ experimentil. 

Neutrino oscillations can occur only if neutrinos are massive particles, if their 
masses are different and if neutrino mixing is realized in nature. In this case, the left- 
handed flavor neutrino fields VaL = e, /i, r) are superpositions of the left-handed 
components UkL {k = l,...,n) of the fields of neutrinos with definite masses m^: 
I'aL = X]fc=i ^ak ^kL , whcrc U IS a Unitary mixing matrix. The general expression for 
the probability of — > transitions in vacuum is 
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where Amlj = ml — m^, L is the distance between the neutrino source and detector 
and p is the neutrino momentum. 

The three experimental indications in favor of neutrino oscillations require the 
existence of three different scales of neutrino mass-squared differences: Am^^j^ ~ 
10-5 eV^ (MSW) or Am^^^ ~ lO-^^eV^ (vacuum oscillations), Am^^^ ~ 5 x IQ-^eV^ 
and AmggL ~ 1 eV^, where Amgg^ is the neutrino mass-squared difference relevant for 
short-baseline (SBL) experiments, whose allowed range is determined by the positive 
result of the LSND experiment. 

Three independent mass-squared differences require at least four massive neu- 
trinos. The number of active light flavor neutrinos is known to be three from the 
measurement of the invisible width of the Z-boson, but there is no experimental up- 
per bound for the number of massive neutrinos (the lower bound is three). In the 
following we consider the simplest possibility of existence of four massive neutrinos. 
In this case, in the flavor basis, besides the three light flavor neutrinos Ue, z/^, Ur that 
contribute to the invisible width of the Z-boson, there is a light sterile neutrino Ug 
that is a S\J{2)l singlet and does not take part in standard weak interactions. 

Two years ago we have showni that among all the possible four-neutrino mass 
spectra only two are compatible with the results of all neutrino oscillation experi- 
ments: 

atm sun sun atm 

(A) mi < 7712 ^ ms < 7714 and (B) m-i < 7712 ^ m-a < 7714 . (2) 
^ ■' >> ^ ■' 

SBL SBL 

In these two schemes the four neutrino masses are divided in two pairs of close 
masses separated by a gap of about 1 eV, which provides the mass-squared difference 
AmgBL — ^^41 = ^4 ~ ""^1 that is relevant for the oscillations observed in the LSND 
experiment. In scheme A Am^^^^ = Am2i = m\ — m\ is relevant for the explanation of 
the atmospheric neutrino anomaly and Am^^^^ = Am43 = 7n\ — m| is relevant for the 
suppression of solar z/e's, whereas in scheme B Am^^j^ = Am43 and Am^^^ = Am|^. 
Let us define the quantities c^, with a = e, /i, r, s, in the schemes A and B as 

(A) c^=Y,\Uc.k\\ (B) c^=Y,Pc.k?- (3) 

k=l,2 fc=3,4 

Physically Cq, quantify the mixing of the flavor neutrino with the two massive 
neutrinos whose Am^ is relevant for the oscillations of atmospheric neutrinos (z/i, z/2 
in scheme A and z/3, z/4 in scheme B). 

The probabihty of transitions (/3 7^ a) and the survival probability of z/^ 

in SBL experiments are given byi 



with the oscillation amplitudes 



2 

Ba;a = 4C„ (1 - Ca) , (5) 

where the index k runs over the values 1,2 or 3,4. The probabilities (^) have the 
same form as the corresponding probabilities in the case of two-neutrino mixing, 
Pt,^^up = sm'^{29) sin2(Am2L/4p) and Pu^^u^ = 1 - sin2(26') sin^ {Am'^ L / ip) , which 
have been used by all experimental collaborations for the analysis of the data in order 
to get information on the parameters sin^(2^) and Am^ {9 and Am^ are, respectively, 
the mixing angle and the mass-squared difference in the case of two-neutrino mixing). 
Therefore, we can use the results of their analyses in order to get information on the 
corresponding parameters Aa-p, Ba-a and AmggL- 

The results of all neutrino oscillation experiments are compatible with the schemes 
A and B only ifl 

Ce < Og and > 1 — a° , (6) 

where 

«° (« = e,/i) (7) 

and B^.^ is the upper bound for the amplitude of ^i^a^^^a oscillations obtained from 
the exclusion plots of SBL reactor and accelerator disappearance experiments. Hence, 
the quantities Og and a° depend on Am^. The exclusion curves obtained in the Bugey 
reactor experiment and in the CDHS and CCFR accelerator experimental imply that 
both a° and a?, are smalli: a° < 4 x 10~^ and < 2 x 10~^ for any value of Am^ in 
the range 0.3 < Am^ < lO^eV^. 

The smallness of Ce in both schemes A and B is a consequence of the solar neutrino 
problemE. It implies that the electron neutrino has a small mixing with the neutri- 
nos whose mass-squared difference is responsible for the oscillations of atmospheric 
neutrinos (z/i, in scheme A and z/3, z/4 in scheme B). Therefore, the transition prob- 
ability of electron neutrinos and antineutrinos into other states in atmospheric and 
long-baseline (LBL) experiments is suppressed. Indeed, it can be shownia that the 
transition probabilities of electron neutrinos and antineutrinos into all other states 
and the probability oi^v^-^Vf. transitions in vacuum are bounded by 

1 - P:^^1\ < a° (2 - a°) (8) 
i^!^^,^J,<minfa°(2-a°),a° + iAO (9) 

where A°.g is the upper bound for the amplitude oi'^v^.—^Vf. transitions measured in 
SBL experiments with accelerator neutrinos. 
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Figure |l] 

The two schemes A and B have identical imphcations for neutrino oscillation ex- 
periments, but very different implications for neutrinoless double-/? decay experiments 
and for tritium /3-decay experiments. Indeed, in scheme A 

\{m)\ < nii and m('^H) ~ m4 , (10) 

whereas in scheme B 

\{m)\ < a'^mi nii and m(^H) <^ m4 , (11) 

where (m) = Ylt=i ^ei effective Majorana mass that determines the matrix 

element of neutrinoless double-/? decay and m(^H) is the neutrino mass measured in 
tritium /?-decay experiments. Therefore, in scheme B \{m)\ and m(^H) are smaller 
than the expected sensitivity of the next generation of neutrinoless double-/? decay and 
tritium /5-decay experiments. The observation of a positive signal in these experiments 
would be an indication in favor of scheme A. 

Summarizing, the results of neutrino oscillation experiments indicate that only 
the two four-neutrino schemes (^) are allowed and the electron neutrino has a very 
small mixing with the two massive neutrinos that are responsible for the oscillations 
of atmospheric neutrinos {ui, V2 in scheme A and 1/3, 1/4 in scheme B). Hence, the 
two schemes have the form shown in Fig.|l], where is associated with the two 
massive neutrinos neutrinos that are responsible for the oscillations of solar neutrinos 
(z/3, 1/4 in scheme A and z^i, z/2 in scheme B), with which it has a large mixing, 
whereas i/^ is associated with the two massive neutrinos neutrinos that are responsible 
for the oscillations of atmospheric neutrinos, with which the muon neutrino has a 
large mixing. The results of neutrino oscillation experiments do not provide yet an 
indication of where Vr and Vg have to be placed in the two schemes represented in 
Fig.|l[ We will show in the following that the standard Big-Bang Nucleosynthesis 
constraint on the number of light neutrinos provide a stringent limit on the mixing 
of the sterile neutrino with the two massive neutrinos that are responsible for the 
oscillations of atmospheric neutrinosi'0. 

It is well known that the observed abundance of primordial light elements is 
predicted with an impressive degree of accuracy by the standard model of Big-Bang 



Nucleosynthesis if the number A^^ of hght neutrinos (with mass much smaller than 
1 MeV) in equilibrium at the neutrino decouplingtemperature {T^cc — 2 MeV for z/g 
and Tdec — 4 MeV for u^, u^) is not far from threeliiHi3. 

The value of is especially crucial for the primordial abundance of ^He. This is 
due to the fact that N^y determines the freeze-out temperature of the weak interaction 
processes + n ^ p + ij^, Ug + n ^ p + e~ and n ^ p + + Ue that maintain protons 
and neutrons in equilibrium, i.e. the temperature at which the rate TwiT) ~ GpT^ 
{Gp is the Fermi constant) of these weak interaction processes becomes smaller than 
the expansion rate of the universe 



(Mp is the Planck mass), where R{T) is the cosmic scale factor and g^, = 5.5 + 
1.75Ni, for me ^ T < m^. If N^, = 3 the primordial mass fraction of ^He, Yp = 
mass density of '^He / total mass density, is Yp ^ 0.24, which agrees very well with 
the observed valuelli Yp^^ = 0.238 ± 0.002. Since Yp is very sensitive to variations of 
Niy, it is clear that the observed value of Yp puts stringent constraints on the possible 
deviation of A^^ from the Standard Model value N^, = 3. 

In the four-neutrino schemes (0) standard BBN gives a constraint on neutrino 
mixing if the upper bound for N,^ is less than four. In this case the mixing of the sterile 
neutrino is severely constrained because otherwise neutrino oscillations would bring 
the sterile neutrino in equilibrium before neutrino decoupling, leading to Ni, = 4. In 
particular, we will show that standard BBN with < 4 implies that Cs is extremely 
small. 

The amount of sterile neutrinos present at nucleosynthesis can be calculated using 
the differential equationlil 



where 71,^^ is the number density of the sterile neutrino relative to the number density 
of an active neutrino in equilibrium and r^,^ are the collision rates of the active 
neutrinos, including elastic and inelastic scattering^, T^^ = 4.0 G^pT^ and T^^ = 
Ti,^ = 0.7 r,^^. The quantities {Puc^-^uJ coii are the probabilities of Ua — > transitions 
averaged over the collision time tcou = ^/^u^- Hence, also n^^ has to be considered as 
a quantity averaged over the collision time. 

Equation ([13|) describes non-resonant and adiabatic resonant neutrino transitions 
if ^osc -C tcou texp- The condition tosc -C tcou means that neutrino oscillations have 
to be fast relative to the collision time. The characteristic expansion time of the 
universe texp is given by tgxp = ^/H where H is the Hubble parameter H = R/ R, 
which is related to the temperature T by if = —T/T ~ 0.7 (T/1 MeV)^ s^^ (this 
value can be obtained from Eq.(|12D with me ^ T < m^ and 3). The relation 
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Ty^l H ~ 1.2 (T/1 MeV)^ shows that for temperatures larger than 2 MeV the colhsion 
time is always much smaller than the expansion time@. 

Since by definition Ny is the effective number of massless neutrino species at Tdec, 
in order to get a constraint on the mixing of sterile neutrinos we need to calculate 
the value of n^^ at Tjec produced by neutrino oscillations. With the initial condition 
n^XTi) = {Ti ^ 100 MeV), the integration of Eq.(|T|) gives§ n^^(Tdec) = 1-6"^ 
with ^ 

^ Yl {Pu^^uJ coiiT u^dT . (14) 

'^dcc a=e,fj.,T 

Imposing the upper bound n,^^(Tdec) < SN = Nu — 3 one obtains the condition 
F < \ \n{l -SN)\. 

For the calculation of F the averaged transition probabilities (Puc^us) coii must be 
evaluated and the effective potentials of neutrinos and antineutrinos due to coherent 
forward scattering in the primordial plasmalll, 



'4 



V, = -Q.02Gfp^ = V, V^,r = ^V and K = 0, (15) 

must be taken into account (in the absence of a lepton asymmetry the effective po- 
tentials of neutrinos and antineutrinos are equal). Here p is the neutrino momentum, 
which we approximate with its temperature average (p) ^ 3.15 T, Gp is the Fermi 
constant, Mw is the mass of the W boson and ^ = cos"^ 9^/(2 + cos^6'p^) ~ 0.28, 
where 6w is the weak mixing angle. The propagation of neutrinos and antineutrinos 
is governed by the effective hamiltonian in the weak basis 



Hw =P + diag 
2p 



m^^, mg, 



[/t + diag V,^V,^V,0 . (16) 



It is convenient to subtract from Hw the constant term p + mi/2p + , which does 
not affect the relative evolution of the neutrino flavor states, in order to get 



H' 



w 



2p 



U diag 



0, Am^i, Am^i, Am^i f/^" + diag {1 - OV,0,0, -^V 



(17) 



From this expression it is clear that the relative evolution of the flavor neutrino 
states depends on the three mass-squared differences and not on the absolute scale 
of the neutrino masses. The effective hamiltonian in the mass basis is given by 



M 



H'm = ^ diag 



0, Am2i, Amg-^, /S.m\^ 



dia£ 



u . 



In the mass basis the mixing has been transferred from the mass term to the potential 
term. In order to calculate the evolution of the neutrino flavors it is necessary to 
parameterize the 4x4 neutrino mixing matrix U . However, since the second and 
third rows and the second and third columns of the diagonal potential matrix in 
Eq.([T8|) are equal to zero, it is clear that the values of the second and third rows of 
U , corresponding to and v-r-, are irrelevant and do not need to be parameterized. 
Furthermore, since Ce is small in both schemes A and B, it does not have any effect 
on neutrino oscillations before BBN and the approximation Cg = is allowed. Hence, 
the 4x4 neutrino mixing matrix in scheme A can be partially parameterized as 



U 











cos^ 



sin 6' \ 



(19) 



y sin if sin x ~ sin (y9 cos x ~ cos (p sin 9 cos >p cos 9 j 



with < < 7r/2. The partial parameterization of the mixing matrix in scheme B 
can be obtained from Eq.([T^) with the exchanges 1^3 and 2 ^ 4 of the columns 
of U . In this way Cg = sin^ in both schemes A and B. The dots in Eq.(|l9|) indicate 
the elements of the mixing matrix belonging to the and u^- rows {U^i and Uri 
with i = 1, ...,4), which do not need to be parameterized. In Eq.(p!9|) we have 
parameterized only the elements of the mixing matrix belonging to the and Ug 
lines {Uei and Ugi with i = 1, . . . , 4) in terms of the three mixing angles 9, Xi ^- If is 
clear that this partial parameterization of the mixing matrix (with the approximation 
Uei = Ue2 = 0) is much easier to manipulate than a complete parameterization, which 
would require the introduction of 6 mixing angles and 3 complex phases. 

Notice that no complex phase is needed for the partial parameterization of the 
mixing matrix in Eq. ([l9|) , because the elements Uei and Usi with i = 1, ... ,4 can 
be chosen real. Indeed, the line Usi with i = 1, ... ,4 and the element Ue4 can be 
chosen real because all observable transition probabilities are invariant under the 



phase transformation [/, 



where Xa and Uj are arbitrary parameters. 
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In order to show that also t/es can be chosen real, we multiply the unitarity relation 
Y2t=i U*f. = with f/g*3 Us3- The imaginary part of the resulting relation gives 



Im[[/e3 U:, U:, Us4 = Im[f/ei U:, U:, U,s] + lm[Ue2 U:, U:, Us3] 



(20) 



In the approximation Uei = Ue2 = the right-hand part of Eq.(^) vanishes. There- 
fore, since t/^s, ?7e4, Usi have been chosen to be real, also f/es must be real. 

Since the mass-squared differences have a hierarchical structure, Am^g ^ Ar/igi ^ 



Am\'^ in scheme A and Am,! Am^ 
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43 



in scheme B, the effective hamilto- 



nian H'j^ can be diagonalized approximately taking into account only one of the three 
Am^'s for different ranges of the temperature T. Then, it can be shown thatilli the 
condition F < | ln(l — 5N)\ gives the bound 



920 



leV^ 



1/2 



ds a/1 - 4 + 33 



Am 



atm 



1/2 



sin^ 2x 
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10-2 eVV VI + cos 2x 



< 



ln(l-5iV)|, 
(21) 

with dg = Cg in scheme A and dg = 1 — Cs in scheme B. 

Both terms in the left-hand side of Eq.(|2l|) are positive and must be small if 
6N < 1. The SBL term, depending on AmggL, is small if either Cg is small or large, 
but the atmospheric term, which depends on Am^^^^, is small only if Cg is small. 
Indeed, if Cg is close to one we have {U^i,Ufj_2) ~ (cos x, sin x) in scheme A and 
{U^s,U^4) ~ (cos X, sin x) in scheme B. This means that, in order to accommodate 
the atmospheric neutrino anomaly, sin^ 2x cannot be small. This is in contradiction 
with the inequality (^) and we conclude that the bound (pT]) implies that Cg is small. 

Since Cg is small only non-resonant transitions of active into sterile neutrinos due 
to AmgBL are possible in scheme A, as illustrated in Fig.^ where we have plotted 
the effective squared masses (obtained from a numerical diagonalization of the hamil- 
tonian (pIBf)) as functions of (z/g does not have resonant transitions into z/^j or Uj. 
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because we have chosen Cg = 0). Hence, the conditions for the vahdity of Eg .([131) are 
satisfied and the SBL term in Eq. (^TD gives the bound 

-1/2 

|ln(l-5Ar)|. (22) 



/\t7?,^ 

c, < 1.1 X 10"^ ' 



leV 

On the other hand, since Cg is small, a resonance occurs in scheme B at the tempera- 
ture Tres = 16(Am|BL/l eV^)i/^|l -2c,|^/^ MeV, as illustrated in Fig.|. The condition 
5iV < 1 implies that this resonance must not be passed adiabatically. In this case the 
conditions for the validity of Eq. ([l3|) are not fulfilled and the SBL term of Eq. (|2l|) 
does not apply. Using an appropriate formula0 for the calculation of the amount 
of sterile neutrinos produced at the resonance through non-adiabatic transitions one 
show0 that the BBN bound on in scheme B is given by 
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leV^ 
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ln(l -5A^)| . 



(23) 



Figure ^ shows the values of the bounds (^) and (p3|) obtained from the LSNDB 
lower bound Am^BL ^ 0.27 eV^ for 0.2 < < 1. One can see that standard BBN 
implies that is extremely small. Therefore, is mainly mixed with the two mas- 
sive neutrinos that contribute to solar neutrino oscillations (z/3 and z/4 in scheme A 
and Vx and in scheme B) and the unitarity of the mixing matrix implies that 
is mainly mixed with the two massive neutrinos that contribute to the oscillations of 
atmospheric neutrinos. Adding this information to the two schemes depicted in Fig.|l] 
we obtain the schemes shown in Fig.^. These schemes have the following testable im- 
plications for solar, atmospheric, long-baseline and short-baseline neutrino oscillation 
experiments: 



The solar neutrino problem is due to Vf. 



oscillations. This prediction will 



be checked by future solar neutrino experiments that can measure the ratio of 
neutral-current and charged-current eventsEl. 



The atmospheric neutrino anomaly is due to u^^ 
will be checked by LBL experiments. 



Ur oscillations. This prediction 



• ^ Ur and Ug oscillations are strongly suppressed in SBL experiments. 

With the approximation Cg — 0, for the amplitude of oscillations we have 

the upper bound A^-r < {(^eY, that is shown in Fig.^ (solid curve) together with a 
recent exclusion curve obtained in the CHORUS experiment (dash-dotted curve) 
and the final sensitivity of the CHORUS and NOMAD experiments (dash-dot- 
dotted curve)Ei. 

If these prediction will be falsified by future experiments it could mean that some 
of the indications given by the results of neutrino oscillations experiments are wrong 
and neither of the two four neutrino schemes A and B is realized in nature, or that 
Big-Bang Nucleosynthesis occurs with a non-standard mechanism^. 

In conclusion, we would like to emphasize that if the analysis presented here is 
correct and one of the two four neutrino schemes depicted in Fig.^ is realized in 
nature, at the zeroth-order in the expansion over the small quantities Ce and Cg the 
4x4 neutrino mixing matrix has an extremely simple structure in which the z/g, Ug 
and z/^, Uj. sectors are decoupled. For example, in scheme A 



U 



cos^ sin^ \ 

cos 7 sin 7 

— sin 7 cos 7 

\ -sin^ cos^ / 



(24) 



where 6 and 7 are, respectively, the two-generation mixing angles relevant in solar and 
atmospheric neutrino oscillations. Therefore, the oscillations of solar and atmospheric 
neutrinos are independent and the two-generation analyses of solar and atmospheric 
neutrino oscillations yield correct information on the mixing of four-neutrinos. 
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